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1 Introduction 



The problem of finding a classical r-matrix structure for the Calogero-Moser (CM) type 
of models aroused some attention a few years ago, cf. The fact that this had 

remained an open problem until relatively recently lies probably in the specific feature 
that for these models the r-matrix turns out to be of dynamical type, i.e. it depends on 
the dynamical variables. Similar features have been found in other integrable many-body 
problems as well, e.g. systems separable in the generalized ellipsoidal coordinates [Q]. The 
difficulty presented by the dynamical aspect of the r-matrix is that the Poisson algebra 
of a model, whose structural constants are given by a dynamical r-matrix is, generally 
speaking, no longer closed, and that there is no closed-form Yang-Baxter equation defin- 
ing the r-matrix. So far, only for one particular example — the spin generalisation of 
the Calogero-Moser model — a proper algebraic setting (Gervais-Neveu-Felder equation) 
is found || which also allows to quantize the model. For other models finding the alge- 
braic interpretation of the dynamical r-matrix and, respectively, solving the quantization 
problem are still open questions. Thus, the use of such dynamical r-matrices so far is 
very restricted. Nonetheless, the existence of such structures is probably significant for 
the integrability of the model — even though one has not been able thus far to deal very 
effectively with these structures — and it is foreseeable that they will play a role in the 
further understanding of the models, maybe even on the quantum level. 

One of the most important integrable many-body systems is the relativistic variant of 
the Calogero-Moser model, the so-called Ruijsenaars- Schneider (RS) model introduced in 
H 0]. Its importance lies in the fact that it can be considered as a g-deformation of the CM 
model and as such the corresponding quantum model is realised in terms of commuting 
difference operators whose eigenfunctions are given in terms of Macdonald polynomials, cf. 
e.g. [H, |9| . On the classical level a dynamical r-matrix was found only very recently in ]l0| 
for the rational and trigonometric (hyperbolic) cases, although a special parameter-case 
was already treated in an earlier paper, ill]. A geometric interpretation was given in a 



recent preprint, cf. [12|. So far, no results have been found for the full elliptic case. That, 
in fact, is the subject of the present paper where we will present the dynamical r-matrix 
structure for the RS model in the generic elliptic case, thus generalising the previous 
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results of [10 1 



2 Ruijsenaars- Schneider model and its discrete-time version 

The equations of motion of the RS model in its generic (elliptic) form read 

q\ = QiQjv(qi - qj) , i=l,...,N , (2.1a) 
where the potential v{x) is given by 

v(x) = ( ^' {X \ ; , (2.1b) 
p(A) - p[x) 

in which p(x) = p(x\u>i,u>2) is the Weierstrass P-function, 2wi t 2 being a pair of periods, and 
A is the (relativistic) deformation parameter. As shown by Ruijsenaars and Schneider in 
H 0], this multi-particle model is integrable, and carries a representation of the Poincare 
algebra in two dimensions. Moreover, a large number of the characteristics of the CM 
model are generalized in a natural way to the relativistic case, such as the existence of a 
Lax pair, a sufficient number of integrals of the motion in involution, and exact solution 
schemes in the special cases of rational and trigonometric/hyperbolic limits. The elliptic 
case has recently been investigated by Krichever and Zabrodin in [|l3| in connection with 
the non-abelian Toda chain. 



In [14| there was constructed an exact time-discretization of the equations (2.1a) given 



by an integrable correspondence of the form 

N , N <?{qi-qk) cr{qe-q k +X) 

u a(q e -q k + X ~ , i=l,...,N. (2.2) 

tJi v{qe - Ik - A) ^ a{qe - q k ) a{qi - q k - X) 

In ( |2.2p the q k denote the particle positions for the time variable equal to n, the tilde being 
a shorthand notation for the discrete-time shift, i.e. for qk(n) = q k we write q k {n + l) = q k , 
and q k {n — 1) = q k . The function cr(x) is the Weierstrass sigma-function, (see Appendix 



for the definition), and A is the parameter of the system as in the continuous case (2.1). 

The initial value problem for eqs. (|2.2|), given initial particle positions {^(0)} and 
{%(!)}, leads to the problem of solving at each iteration step a coupled system of N 
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algebraic equations for N unknowns, and it was shown in [|14[ that in fact it is an integrable 
symplectic correspondence (for a definition, see e.g. [jl5[) with respect to the standard 
symplectic form Q = J2 k dp k A dq k . This implies that any branch of the correspondence 
given by eqs. (|2.2| ) defines a canonical transformation with respect to the standard Poisson 
brackets given by 

{Pk,Qe} = he , {Pk,Pe} = {Qk,qe} = ■ (2.3) 

Here 

log \a (q e - q k )\ + log \a(q e - q k +\)\) . (2.4) 



N ( 
k=l \ 



The discrete equations of motion ( p.2j ) arise from a discrete Lax pair of the form 

N 

L K = ~ 1j + X ) e ij > ( 2 ' 5a ) 

M * = 12 h 2 i§M - Qj + X)eij , (2.5b) 
using the discrete Lax equation 

l K M K = M K L K . (2.6) 



Notice here that in (|2.q ) we use a different gauge from the symmetric one used in [14]. In 
eqs. ( |2.5| ) the variable k is an additional spectral parameter, and the matrices e^- are the 
standard elementary matrices whose entries are given by {e%j)kl = ^ik^ji- The function $ K 
is called the Baker function and is defined as 

..(,) S + ^ , (2.7) 

which obeys a number of functional relations listed in the Appendix. The auxiliary vari- 
ables ft| can be expressed in terms of the canonical variables, we obtain 

^ <r(q t -q k ) 

In terms of these variables we have the following Poisson brackets: 
Uk, qe} = , {log h 2 k , q e } = 5 M , 

{log h\ , log /i, 2 } = C(q k -q e + X)+ C(q k - q £ - X) - 2((q k - q t ) , k^i. (2.9) 

It is easy to see that in terms of the canonical variables pi and q$, the Lax matrix L K in 
( |2.5a| ) is exactly the same as the one of the continuous RS model, cf. [16]. In fact, taking 



the continuum limit on the discrete-time part of the Lax pair (2.5), namely the matrix M R 
( |2.5b| ), we obtain a Lax pair for the continuous RS model given by equations ( p. la ). Since 
the L K -matrix for the discrete and continuous models is the same, the proof of involutivity 
of the invariants (integrals) Ii = trL~. is the same in both cases, and sufficient to assess 
the Liouville integrability both discrete as well as continuous. The proof can be found 
in the original paper of Ruijsenaars, 0, but is rather involved. Having at one's disposal 
an r-matrix structure would make the involutivity manifest. So far such an r-matrix has 
not been found in the full elliptic case. We will proceed now to establish this r-matrix 
structure. 



3 Classical r-matrix structure 

As was noted recently by Suris, cf. [12], the main difference between the r-matrix struc- 
tures of the relativistic and non-relativistic CM models resides in the fact that the latter 
is given in terms of a linear Lie-Poisson structure/bracket, whereas the former is given in 
terms of a quadratic bracket, cf. also pi] . The Poisson structure for the RS model will 
thus be given in the following quadratic r-matrix form (cf. []l7| , [Igfl ) 

{ L K f L K , } = L K <g) L R ir~ K , - r+ K ,L K <g> L K , 

+ {L K ® 1) 8 + K , (1 ® L K .) - (1 <g> L K i) s~ K , (L K <g> 1) . (3.1) 

The following symmetry conditions must hold for the r-matrices: r^ K , and s^ K , 

Pr% K >P = -r±,, K , PstyP = V, K , (3.2a) 

where P is the permutation matrix in the tensor product of two matrices, i.e. P(A®B)P = 
B (g> A, as well as the condition 

r ty ~ s k,k' = r K,K,' ~ s k,k' ' (3.2b) 

in order that the quadratic Poisson algebra generates Hamiltonian flows for the invariants 
of the model, cf. ||, [l{J. The condition was also formulated in |l9|] in order to 



allow for a quadratic algebra on the lattice in terms of a local Lax representation to be 
integrated to a quadratic algebra in terms of the monodromy matrix. 



The choice of a gauge for the Lax matrices seems to be quite important in that it 
influences to a great extent the complexity of the associated r-matrix. The Lax matrix 



( 2.5a ) has the nice property that it yields a remarkably simple r-matrix structure even in 
the elliptic case. In fact, the r-matrices we found are of the form 



S Zk> = S K +U 



where [] 



r K,n' = r° >K > + ^ C(« - K')eu <g> en + ^ Cili ~ Qj)e-u ® Zjj > 

i i^j 

r K, K ' = Y - Qj) e ij ® eji , 



« = Efe l5 A^)..e 



h.l 



Y C(Qj - Qi ± X)eu <g) e# . 



(3.3a) 

(3.3b) 
(3.3c) 
(3.3d) 

(3.4a) 
(3.4b) 
(3.4c) 
(3.4d) 



The matrix elements in ( |3.4c| ) can be calculated explicitly using the formula ( A.13j ) for the 
inverse of the elliptic Cauchy matrix as well as making diligent use of the elliptic Lagrange 
interpolation formulae ( A. 10 ) and ( A.ll| ), and this yields the following expression 



LZ 1 d\L, 



C(« + NX) - C(A) + J2 (C(ft " Ik - A) - C(« - Qk)) 

k^i 



+(i - <%) 



N 

n 



cr(% - Qk ~ X) 



" = i cr(Qi ~ Qk) 



k 
k^i 



v(Qj ~ Qk) 



N 

33 cr(<?j - 9fc - A) 



$ K +JVA(<?i - Qj) • (3.5) 



The proof of the r-matrix structure ( |3.1| ) together with ( p.3| ) and ( |3.4[ ) is by direct compu- 
tation starting from the explicit form of the L-matrix Q2.5a|) and the Poisson brackets ( p.S| ) 



1 In (3.4c) we mean by the matrix d\L K the following 

JV 

d\L K = hi$ K (q.i - q 3 + A) [((k + q t - q 3 - + A) - - + A)] ey 

i,j=l 



i.e. we differentiate only with respect to the explicit dependence on the parameter A. 
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and making use of a number of elliptic relations which are listed in the Appendix. We 
will not give any details, but just restrict ourselves to giving a few intermediate relations, 
which can be established using the formulas from the Appendix, namely: 



(L K <g) L K >) s K = Y Y h i h i'^K(qi - Qj + A) - qy + A) ey ® e V y 



i j v y 



x %' [C(« + Qi- qj + A) - C(q% - qj + A)] , 

(L K <g> 1) s K (1 (g> L K /) = X! X! hihi,® K (qi - qj + A) - + A) <g> e^/ 

U i'j' 



(3.6a) 



x [c(k + % - gj + A) - C(g« - qj + A)] , 



(3.6b) 



as well as 



r K, K > > L « ® = Y Y hihl® K (qi - qj + A) - g^ + A) e„ e^/ 

x { (1 - 5jj,)(l - [C(g» - &/) + C(gi' - qj + A) + C(g, - qy) - C(% - ?j' + A)] 
+ (1 - 6u') [C(qi - qv) - C(« + g« - gj + A) + + g^ - gj' + A) + £(« - «')] 

- (5jj/) [C(gj - gj') + C(« + gi - g^ + A) - C(«' + gi' - gj' + A) - £(« - «')] } . 

(3.6c) 

Remark here that our r-matrices do not depend on momenta, like in the non-relativistic 
case ||, which was the motivation for the choice of the gauge of L K . 

As a direct application of the r-matrix structure let us calculate the (continuous) time 
part of the Lax representation. It is obtained from the following formula: 



(tr ® id) (L K ®l)(r+ 



(3.7) 



where 



a^)hf + y ^aqi - qj) - Y - % - A ) 

j^i 3 



(3.8) 



which, together with ( 2.5a| ), leads to the Lax representation found in [16] for the continuous 
RS model (up to a gauge transformation!). Thus, the continuous equations of motion ( 2.1a| ) 
corresponding to the Hamiltonian tr L K follow from the Lax equation 



(3.9) 
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Remarks: 



• The non-relativistic limit is obtained by letting A — > while scaling the momenta 
Pi := Xpi and making the canonical transformation p$ := p, L + J2k^i C(<2t ~~ Ik) such 
that h% — ► 1 + Apj + 0(A 2 ) in ( p.8[ ). The r-matrix structure is linear in that limit 
since the L-matrix behaves as 

L K -> A" 1 + C(k) + ^Pie^ + - qj)eij + 0(A) , 

i i^j 

whereas the matrices r„ „/ , . enter in the following combination 

rt, K > + s- K , - r«$ + 0(A), (3.10) 
in which the non-relativistic r-matrix is given by 

+ ^ - grj)ejj ® ey , (3.11) 

thus recovering the result of [|| in the leading terms. 



It does not seem easy to compare our result with that of Avan and Rollet in [10] 
because of the not so transparent nature of that result (they wrote their r-matrix 
structure in a linear form hiding the quadratic nature of the RS model). In the 
trigonometric limit taking the special parameter value e 2A = — 1 the structure we 
found should reduce to the one given earlier by Babelon and Bernard in [11 1. 



In []13] the author seems to suggest that the r-matrix for the relativistic model and 
the non-relativistic model is the same in the rational and trigonometric/hyperbolic 
limits, which is demonstrated by an independence of the r-matrix objects on the 
relativistic parameter A. This, however, seems to be no longer true in the case of 
the elliptic model, and there is no apparent way in which one can get rid of the 
A-dependence in the explicit formulas (|3.3|) and 



• We do not write down any Yang-Baxter type relations between the r-matrices given 
in (|3.3|), because as a consequence of the dynamical nature the Yang-Baxter algebra 



does not seem to be closed, i.e. the Yang-Baxter 2-cocycle consists of terms which 
contain Poisson brackets with the L-matrix itself. It is an interesting open problem to 
see whether one can close the algebra on any level, leading to a possible truncation 
of some higher-order Yang-Baxter cocycle. So far, no results along this direction 
exist. 



4 Conclusions and outlook 

In this paper we have presented the classical r-matrix structure of the full elliptic Ruij- 
senaars-Schneider model. Since the model is the most general among the Calogero-Moser 
type models for the sl n Lie algebra, our result is in a sense conclusive. Still, a number of 
questions wait to be answered. Since the dynamical nature of the r-matrices implies that 
the corresponding Yang-Baxter algebra is not closed, it is not yet clear how to use it for 
quantization. 



Concerning the quantization problem, the recent result of Hasegawa [20] should be 
mentioned, who has found an interesting connection between the quantum L-operator 
associated with Belavin's .R-matrix and the quantum integrals of the Ruijsenaars' model. 
This somehow implies that on the classical level there should exist a gauge transformation 
involving the dynamical variables between an elliptic r-matrix of Belavin type and the one 
we have constructed in the present paper. 

Another result concerns the separation of variables approach leading to the explicit 
integral representations for the Macdonald polynomials associated with the trigonometric 
RS model, cf. |2l|] . So far the only result for the elliptic case is the separation of variables 



for the 3-particle nonrelativistic Calogero-Moser model [22|. The elliptic r-matrix could 



presumably help in constructing a separation of variables in the general case. 

One more possible application of the results of this paper could lie in the time dis- 



cretization of the RS model constructed in [14|. One feature of the proposed time dis- 



cretization is that these discrete models share the time-independent part of the Lax pair 
with the corresponding continuous models and, consequently, the invariants take the same 
form in both cases. Thus, the proof of the Liouville integrability (or the involutivity of 
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the invariants) is exactly the same as for the continuous model. For the discrete models 
it seems that the most prominent role is played by the M-matrix. In fact, in similar sys- 
tems related to integrable lattices it was found that there exists an extended Yang-Baxter 
structure which incorporates the M-matrix as well as the L-matrix in the Yang-Baxter 
algebra. Now that a classical r-matrix structure is available for the generic RS system, 
it would be interesting to search for similar extended YB structures for these many-body 
systems as well. This would possibly yield some new insights into the problem of perform- 
ing the i?-matrix quantisation for dynamical r-matrices as the one we have obtained in 
the present paper. 
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Appendix: Formulas for elliptic functions 



Here, we collect some useful formulas for elliptic functions, see also the standard textbooks 



e.g. [23]. The Weierstrass sigma-function is defined by 



<r(x) = x Yl ( 1 ) ex P 

(W(o,o) V UmJ 



- + \{— 



(A.l) 



with Uf-i = 2ku% + 2£ui2 and 2wi ) 2 being a fixed pair of the primitive periods. The relations 
between the Weierstrass elliptic functions are given by 



C(x) = — — , p(x) 
a{x) 



<(*) , 



(A.2) 



where cr(x) and C( x ) are °dd functions and p(x) is an even function of its argument. We 
recall also that the <j{x) is an entire function, and C( x ) is a meromorphic function having 
simple poles at wy, both being quasi-periodic, obeying 



C(x + 2ui 2) = ((%) + 2r?i,2 , cf{x + 2oj\ 2 ) = —a(x)e 



2771,2(^+^1,2) 
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in which 771,2 satisfy r\\i02 — r]2^i = y > whereas p(x) is doubly periodic. From an algebraic 
point of view, the most important property of these elliptic functions is the existence of a 
number of functional relations, the most fundamental being 

ft \ 1 fta\ 1 tt \ ft L«j ^ (r(a + (3)a(P + j)a(-f + a) 

cr(a)cj(/5)cr(7)<T(a + P + 7) 

which can also be cast into the following form 

$ K (x)$ K (y) = <5> K {x + y) [C(as) + Cfr) + CO) - C(« + a? + 2/)] > (A.4) 



The following three-term relation for cr(x) is a consequence of (A. 3) 



a{x + y)a{x — y)a(a + 6)cr(a — 6) = cr(ic + a)a(x — a)a(y + 6)cr(y — b) 

— a(x + b)a(x — b)a(y + a)a(y — a) , (A. 5) 

and this equation can be cast into the following convenient form 

3> K (x)$ x (y) = $k{x - y)<5> K +\{y) + <s> K+ x(x)$ x (y - x) , (A.6) 

which is obtained from the elliptic analogue of the partial fraction expansion, i.e. eq. 

There are few additional important identities that are used in the proof of the r-matrix 
structure, the main one is given by 

$ K - K >(a - + &)*«/(a + y) - $ K - K >(x - y)$ K {y + a)$ K ,(x + b) = 

= <& K (x + a)$ K ,(y + b) [C(a - b) + C(x + b) - C(x - y) - ((y + a)] , (A.7) 

which can be derived from ( |A.6| ) together with ( |A.4| ), and 

®k-k>(x - y)^ K {y + a)$> K r(x + a) = 

= $ K (x + a)$ K ,{y + a) [((x - y) - £(« + x + a) + ((k' + y + a) + C(k - «/)] . (A.8) 

It is eqs. ( |A.7| ) and ( |A.8[) that are used in the derivation of ( |3,6c| ) which forms the main 
step in the computation of the r-matrix. 

In there was used an elliptic version of the Lagrange interpolation formula, which 
was derived on the basis of an elliptic version of the Cauchy identity. We can write the 
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elliptic Cauchy identity in the following elegant form: 

J[k<i a ( x k ~ x e )a(y e - y k ) 



det (^(xi -%•)) = *«(E)<7(E)- 



Y\ki a ( x k - ye) 



where E = ^(x* - y,-) . 

(A.9) 



An elliptic form of the Lagrange interpolation formula is obtained by expanding 
along one of its rows or columns. Thus, we obtain 



N 

n 



N 



and 

N 

n 



N 



U.f=i <r(Vi ~ Vj) 



when E = Y,(xi-Vi)¥>0, (A- 10) 



i=l 



UjLi<r(Vi-Bj) 



E[c(^-b)-c(^-k)] 



JV 



i=l 



when ^2( yi - Xi ) = 0. 



i=l 



(A.11) 



(here x denotes one of the zeros Xj). Note that in this case the left hand side is a mero- 
morphic function on the elliptic curve as a consequence of Abel's theorem. It can be easily 



verified that eq. (A. 11) is independent of the choice of x as a consequence of the relation 



J2— 

i=i U*=i <r{Vi - Vj) 



when ^2(yi - Xi) = , 



(A.12) 



i=l 



(cf. e.g. [23], p. 451), which follows from eq. ( [A. 10 ) in the limit E — > 0. 

Finally, we give the expression for the inverse of the elliptic Cauchy matrix, namely 

P(vi)Q{xj) 



($ K (x. -y.)) 1 .. = $ K+S (y; -Xj) 



(with E as before), in terms of the elliptic polynomials 



Qi(yi)Pi{xj) 



(A.13) 



N N 

*;=! fc=i 



and 



(A.14) 



k^j k^i 

Equation (|A.13|) can be derived using ( |A.10[) and ( |A.11[ ), and it is used to derive equation 
( |3.5|) in the main text. 
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